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Abstract. We study the dynamics of a class of two dimensional stochastic 
processes, depending on two parameters, which may be interpreted as two different 
temperatures, respectively associated to interfacial and to bulk noise. Special lines 
in the plane of parameters correspond to the Ising model, voter model and majority 
vote model. The dynamics of this class of models may be described formally in terms 
of reaction diffusion processes for a set of coalescing, annihilating, and branching 
random walkers. We use the freedom allowed by the space of parameters to measure, 
by numerical simulations, the persistence probability of a generic model in the low 
temperature phase, where the system coarsens. This probability is found to decay 
at large times as a power law with a seemingly constant exponent 9 0.22. We also 
discuss the connection between persistence and the nature of the interfaces between 
domains. 
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1. Introduction 

As is well known, an Ising system quenched from high temperature to low temperature 
exhibits coarsening in its temporal evolution Jj], This property holds in any 
dimension. The voter model ||, defined as a purely dynamical system, is identical 
to the Ising model with Glauber Q or heat bath dynamics in one dimension but its 
behaviour progressively departs from that of the latter when dimension increases. For 
instance the two dimensional voter model also exhibits properties similar to coarsening 
|(|, though the way the system coarsens differs, in some respect, from that of the 
two dimensional Ising model. Measuring the persistence probability of the two models 
demonstrates more strikingly the difference in behaviour between them. While the 
persistence probability for the two dimensional Ising model at zero temperature decays 
as t~ e with 8 w 0.22 |?|, ||, it behaves as exp[— A In 2 t] for the two dimensional voter 
model 0,[TT). 

In an attempt to understand the differences in behaviour between the dynamics 
of the two dimensional Ising and voter models, we were naturally lead to introduce a 
whole class of models interpolating between them. The idea stems from the observation 
that the rules for updating a spin in the voter model resemble those used for the Ising 
model, in the Glauber or heat bath algorithms. This class of models is defined in 
any dimension d, but the present study is restricted to two dimensions, where the 
models depend continuously on two parameters, which may be interpreted as two 
temperatures. 

The aim of the present work is to use the freedom allowed by the space of 
parameters in order to investigate the mechanisms by which coarsening and persistence 
continuously change when going from the Ising model to the voter model. We will 
see that physical insight is provided by interpreting the two parameters defining the 
dynamics as two different temperatures, respectively associated to interfacial and 
to bulk noise. We shall also investigate whether the persistence exponent 9, which 
seems constant w 0.22 for the two dimensional Ising model when temperature varies 
|L2| , |l3| , is a universal exponent for the whole class of models, in the low 

temperature region of the plane of parameters. 

In section 2 we first focus our interest on the definition of this class of models. 
We then give a qualitative description of its dynamics, showing the existence of a 
critical line between a low and a high temperature region in the space of parameters 
defining the models (section 3). Section 4 is devoted to the study of persistence in the 
low temperature region. We first measure the fraction of spins which never flipped 
up to time t, along the special line joining the voter model to the zero temperature 
Ising model. We then discuss the question of persistence at finite temperature (section 
4). Finally an appendix is devoted to the dual description of the models in terms of 
reaction diffusion processes. 

After this work was completed, we discovered that this class of models had been 
previously introduced in [ [l6| , which also contains a determination of the critical line 
by a finite size scaling analysis in the stationary state. For completeness, we kept the 
original wording of sections 2 and 3, adding references to this work where appropriate. 

2. Definition of the class of models 

Let us consider a two dimensional lattice of spins <ii = ±1, evolving with the following 
dynamical rule. At each evolution step, the spin to be updated flips with the heat 
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bath rule: the probability that the spin <7j takes the value +1 is P(<7j = 1) = p(hi), 
where the local field hi is the sum over neighbouring sites ^ • aj and 

p(h) = ^(l+t<mh[p(h)h]). (1) 

The functions p{h) and are defined over integral values of h. For a square lattice, 
h takes the values 4, 2, 0, —2, —4. We require that p(—h) = 1 — p(h), in order to keep 
the up down symmetry, hence 0(—h) — /3(h). Note that this fixes p(0) = 1/2. The 
dynamics therefore depends on two parameters 

Pi=p(2), P2=P(4), (2) 

or equivalently on two effective temperatures 

Tl = Wr T2 = wr (3) 

Defining the coordinate system 

2 2 

ti = tanh — , t 2 = tanh — (4) 

J-i J-2 

with < t\,t2 < 1, yields 



with 1/2 <pi,p2 < 1- 

Hereafter we shall interpret Ti and T 2 as two temperatures, respectively associated 
to interfacial noise, and to bulk noise. This should be understood in the following 
sense. Consider the initial configuration where the system is divided by a flat interface 
into two halves, one half with all spins a = 1, and the other one with all spins a = — 1. 
Then, if pi = p\ = 1, i.e. Ti = Tfe = 0, this configuration will not evolve in time, 
neither in the bulk, nor on the interface, since all spins are surrounded by at least 
three spins of the same value. However, if T 2 = 0, then as soon as Ti > 0, spins at the 
interface will flip, while those in the bulk will not. Conversely, if Ti = 0, then spins 
in the bulk will flip if T 2 > 0, while those on the interface will not. (They will later 
do so because of the noise coming from the bulk.) Note that a configuration where 
the system is divided into two halves by a curved interface will always evolve, even if 
P\ = P2 = 1, since p(0) = \. 

Each point in the parameter plane (pi,P2), or alternatively in the temperature 
plane ^2), corresponds to a particular model. The class of models thus defined 
comprises as special cases the Ising model, the voter and antivoter models, as well as 
the majority vote model, the description of which follows. 

The Ising model with heat bath dynamics corresponds to choosing 

0(h) = const. = 0, V/i, (6) 
where is the usual inverse temperature. Hence 

Pi = ^(1 + tanh 20) P2 = i(l + tanh 4/3). (7) 
This corresponds to the line 

^= 1-2^+2^ ° r h=t » (8) 
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when P varies. For example, at zero temperature, one has p\ = P2 = 1, hence 
T\ = T2 = 0. The dynamics is therefore only driven by the curvature of the interfaces 
between domains of equal values of the spin ||. While the Ising model possesses a 
well defined energy for which the heath bath dynamics satisfies detailed balance, none 
of the models outside the Ising line shares these properties [|l6| , hence these models do 
not possess a usual equilibrium description, and their definition is purely dynamical. 

The dynamics of the voter model is defined as follows At each time step, the 
spin to be updated is aligned with one of its neighbours, chosen at random. Therefore 



P{h) = \{l + \l (9) 



3 1 

Pi = 4, P2 = 1 or t t = -, t 2 = 1. (10) 

This definition corresponds to a model with no bulk noise (j>2 = 1, or T2 = 0). The 
noisy voter model is defined as a simple generalization of (||). The spin to be updated 
is now aligned with one of its neighbours, chosen at random, with a probability 7 
j| [| . In other terms 

P {h) = hi + lh). (11) 



Hence 
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Pi = 7^(1 + 7}) P2 = \{l+l). (12) 



This corresponds to the line 



p 2 = 2pi - - or ti = ^— -2 , (13) 



1 , h 
— or ii = ., 

2 1 + t z 2 

when 7 varies from to 1. One may also extend the definition ([ll]) to —1 < 7 < 0, by 
allowing negative values of the coordinates t\ and £2 , or equivalently letting p\ and P2 
be less than 1/2. The model thus defined is known as the antivoter model 0. 

Finally, for the majority vote model || [l7]], spins are aligned with the local field 
(i.e. with the majority of neighbours) with some given probability. More precisely, if 

p(h) = ~ (1 + 5 sign/!) (0<5<1), (14) 
and p(0) = \. Hence 

Pi=P2 = l(l + S), (15) 
therefore the model corresponds to the line 

pi=P2, or h = ^ 2 , (16) 

1 + l 2 



i.e. T2 = 2Ti, when S varies from to 1. 
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Figure 1. Phase diagram. Broken lines correspond to the noisy voter model (V), 
the Ising model (I) and the majority vote model (M). The low temperature phase is 
located in the upper right corner, above the transition line (full line). 



Figure [I] shows, in the (ti^) or (pi,P2) planes, the lines corresponding to the 
Ising, voter and majority vote models. Three other lines deserve attention. Firstly, 
the p 2 = 1 line corresponds to models with no bulk noise (T 2 = 0), hence the dynamics 
is only driven by interfacial noise, defined above. Secondly, the p\ — 1 line corresponds 
to models with no interfacial noise (T\ = 0), hence the dynamics is only driven by bulk 
noise. (In both cases the effect due to the curvature of the interfaces is always present, 
as mentioned above.) For these last models, the local spin aligns in the direction of 
the majority of its neighbours with probability one, if the local field h = 2, i.e. if there 
is no consensus amongst the neighbours. If there is consensus amongst them, i.e. if 
h = 4, the local spin aligns with its neighbours with a probability p2 < 1. Finally the 
transition line between the low and high temperature regions is discussed below. 

See reference |D| for similar definitions. The introduction of the temperatures T\ 
and T2, together with their interpretations as given above, is original to the present 
work. A dual description of the dynamics of the models in terms of reaction diffusion 
processes is given in the appendix. 

3. Phase ordering and the critical line 

In this section we investigate the dynamics of phase ordering for the two parameter 
class of models presented here. For pi and pi given, we let the system evolve, starting 
from a random initial configuration. This study demonstrates the existence of a critical 
line, in the (pi, P2) or (ti, £2) planes, between a low temperature region where clusters 
(or domains) grow indefinitely, and a high temperature region where one only observes 
fluctuations at a finite scale. It also provides a visual illustration of the role of the 
two temperatures T\ and T2, respectively associated to interfacial and to bulk noise. 

Two methods for updating the spins are at our disposal. With sequential updating, 
time is continuous and each spin evolve independently of others, at times chosen 
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randomly according to a Poissonian law. The normalization is chosen such that each 
spin evolves once in a mean unit of time. Therefore, in a simulation of a lattice of N 
spins, one has to choose randomly a spin and update it, N times per time unit. 

Instead of using sequential updating, we adopt a slightly different procedure for 
our simulations, namely a parallel updating of the spins. The lattice is divided into 
two odd and even sublattices, and, during a unit of time, the two sublattices are visited 
in turn, with a systematic update of their spins. This division is needed in order to 
avoid undesirable effects due to the simultaneous update of neighbouring sites. This 
process allows a functional parallclization on the computer. Spins are represented as 
single bits, and it is possible to arrange the spins into computer words, such that 
all spins of a word can be updated simultaneously by global logical operations. This 
method is numerically efficient, allowing for long simulations on large lattices. 

In figures 2-7, we have selected snapshots of the same fragment 256 x 256 of a 
512 x 512 lattice with periodic boundary conditions, for three times t = 8, 64 and 512, 
and six different points in the parameter space of the models. 

Figures [| depict, as reference views, the well known appearance of domain 
growth for the Ising model, respectively at zero temperature, and at T — 0.97 T c . 
The Ising critical point corresponds to t\ = t2 = 1/V2- For T < T c clusters grow 
as while at T c , they only grow as i 1 / 2 , where the dynamical exponent z rs 2.17. 
In the high temperature phase, clusters grow until their sizes reach the equilibrium 
correlation length. 




Figure 2. Snapshots of the evolution for the Ising model at zero temperature 
(pi =P2 = 1). 

This distinction between a low temperature and a high temperature phase can be 
generalized to all models of the two dimensional parameter plane. The transition line 
between the low and high temperature regions, depicted in figure [|, was located by a 
systematic investigation of the parameter plane. 

For instance, figures ||, [|, an d H correspond to an exploration of the line P2 = 1, 
with pi decreasing away from the zero temperature Ising case {p\ = pi = 1). On 
this line no thermal fluctuations occur inside the clusters, so their interiors remain 
forever uniformly black or white. When p\ decreases, the boundaries between domains 
take progressively a fractal appearance, in accord with the fact that interfacial noise 
increases along the line. This is particularly visible on figure ^, which illustrates the 
case of a model close to the voter model, in the low temperature region {jp\ — 0.76, 
P2 = 1), and on figure ||, for the voter model itself (pi = 0.75, P2 — 1). Again 
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Figure 3. Snapshots of the evolution for the Ising model at T = 0.97 T c (pi = 0.8602, 
p 2 = 0.9743). 



one observes growing clusters, with a characteristic size behaving as \ft. A similar 
behaviour is observed for all models such that 0.75 <p\ < l,P2 = 1. 




Figure 4. Snapshots of the evolution for a model close to the voter model, in the 
low temperature region (pi = 0.76, p2 = 1). 

Figure || illustrates the case of a model located in the high temperature region, 
near the voter model (pi = 0.72, p2 = 1). In parallel to what is observed for the 
Ising case above T c , clusters grow until their sizes reach the correlation length. The 
transition to the low temperature region along the line p2 = 1 is therefore only induced 
by the interfacial noise. 

We also display in figure [?] the evolution observed for the majority vote model with 
(pi — P2 — 0.93). The transition takes place for p\ = P2 ~ 0.923, in agreement with 
the value (0.925) given in 0. 

Finally exploring the line pi — 1 , we observe a transition between the low and high 
temperature regimes for p 2 ~ 0.849, in agreement with the value (0.854) given in |l6| ]. 
Figu re p| depicts snapshots of the evolution for the model with (pi — 1, p2 — 0.84). 
See Jl6[ for a determination of the critical line by a finite size scaling analysis in the 
stationary state. 
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Figure 7. Snapshots of the evolution for the majority vote model (pi = p2 = 0.93). 



Figure 8. Snapshots of the evolution for the model with (pi = 1, pi = 0.84). 



4. Persistence 



In this section we address the question of persistence in the low temperature region 
of the (pi,p2) parameter plane depicted in figure |l|. We first present numerical 
measurements of the fraction R(t) of spins which never flipped up to time t, for 
models with no bulk noise, i.e. for = or p2 — 1. The case of persistence at finite 



temperature, i.e. in the present context for T 2 > 0, will be discussed in section 4.2 



4--1. The p2 = 1 line 

We measured persistence in the low temperature region [p\ > 0.75) of the line p2 = 1. 
Figure |9| depicts the results of simulations for times up to 5000, on a lattice 3072 x 3072. 
The fraction R(t) of spins which never flipped up to time t behaves, in the early stage, 
as 

]nR(t) A\n 2 t + B\nt + C + ---, (17) 

then one observes a crossover toward the power law 

R(t) ~ r\ (18) 

with an exponent 6 w 0.22 which seems independent of p\. This behaviour therefore 
interpolates between the two 'pure' cases, namely the zero temperature Ising model 
(pi = pi = 1) where ( [i~8"| ) is observed with the same exponent 9 0.22 and 
the voter model, for which @ holds |[o[ 0. 

The origin of these two 'pure' behaviours can be simply traced back to the two 
driving forces of the dynamics in absence of bulk noise, i.e. respectively to interfacial 
noise for (|i"7|), and to the curvature of interfaces for (|l8|). Indeed, consider first the case 
of the Ising model at zero temperature (Ti = T% = 0) , for which the dynamics is only 
driven by the curvature of interfaces. As a consequence, the only flipping mechanism 
for a spin is the slow sweeping mode of the interfaces, hence persistent spins deeply 
buried inside large clusters disappear slowly, according to ( |l8| ) . Consider now the case 
of the voter model (p\ = 0.75, P2 = 1), starting from a configuration with a planar 
interface separating two regions of opposite spins. This interface thickens gradually, 
because of interfacial noise, taking ever more a fractal appearance. At large times it 
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Figure 9. Log-log plot of the persistence probability R(t) on the line p2 = 1, for 
pi = 1, 0.9, 0.8, 0.77, 0.75, from top to bottom. The inset depicts In R(t) for the voter 
model (pi = 0.75) against In 2 t. The system size is L = 3072. 

becomes difficult to speak any longer of an interface, the thickness of which should be 
considered as infinite. The decay of persistent spins is much faster and follows ( |l7| ) . 

This analysis allows the following interpretation of the results given above for 
persistence of a generic model along the P2 = 1 line, and p\ > 0.75. In the early stages 
of coarsening, domains begin to grow and develop their interfaces. As a consequence, 
the fraction of persistent spins in the interface region is expected to decay more rapidly 
than if they were deep inside a cluster, following (|l7|) since this is the dominant process. 
The duration of this early stage is related to the thickness of the interfaces, which 
vanishes for the Ising model where smooth interfaces are observed, and diverges for 
the voter model where a fractal appearance of the interfaces is observed. After this 
first stage, the fraction of persistent spins will decay according to (|l8|). Therefore, 
along the line P2 = 1 (Tb = 0), the resulting law for the decay of R(t) will be a 
weighted sum of ( [l7|) and (|lli|). The amplitude obtained by fitting the asymptotic 
algebraic tail ( fTcf ) is naturally interpreted as the fraction of persistent spins which are 
deep inside the clusters. This fraction vanishes when one approaches the voter model, 
with an apparent power law {jp\ — 0.75) Q , with a ~ 2.3. This result is consistent with 
the fact that the voter model is critical, according to the criterion of section ^. 

Note that, for any model such that p\ > 0.75, the prescription of suppressing the 
bulk noise by forbidding flips of a spin a = 1 (resp. a = — 1) surrounded by four 
spins with the same value 1 (resp. —1), leads to an algebraic decay of the persistence 
probability. This corresponds to a projection of the model onto the line p\ = 1. 
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4-2. Persistence at finite temperature 

For a two dimensional Ising system in presence of thermal noise (i.e. in the present 
context, for p2 < 1), the number of spins which did not flip up to time t decays 
exponentially. A more satisfactory quantity to investigate is the fraction R(t) of 
space which remained in the same phase, coming back to the original definition of 
persistence ]l8| , [l9| , where 'remaining in the same phase' means remaining 'dry', i.e. 
unswept by an interface. The prescription given in []l2| to determine whether a given 
spin 'remained in the same phase', leads again to algebraic decay of R(t), the fraction 
of 'persistent' spins, when T < T c . Further works |D| [l4|, [ll| seem to indicate that 
the persistence exponent 9 for the two dimensional Ising model is constant in the low 
temperature phase, and approximately equal to 0.22. 




Figure 10. Log-log plot of the persistence probability R(t). Full lines, from top 
to bottom: (pi = l,p 2 = 1): Ising at T = ; (pi = 0.95, P2 = 0.997): Ising 
at T = 0.6 T c ; (pi = l,p 2 = 0.95); (pi = 0.95, p 2 = 0.95): Majority vote; 
(pi = 0.9, p 2 = 1); (pi = 0.8,p 2 = 0.988): Ising at T = 0.8T C ; (pi = l,p 2 = 0.9); 
(pi = 0.876, p 2 = 0.980): Ising at T = 0.9 T c ; Dashed lines, from top to bottom: 
(pi = 0.92, p 2 = 0.92): Majority vote; (pi = 0.854, p 2 = 0.971): Ising at T = T c ; 
(pi = 1,P2 = 0.84); (pi = 0.75, p2 = 1): Voter. The system size is L = 1536. 

Coming back to the class of models under investigation, one may therefore wonder 
what is the value of the persistence exponent outside the Ising line, in the low 
temperature region. Here we use the definition of |12]| , adopting the extensions given 
in []l4| , in order to measure persistence in the low temperature region of the parameter 
plane. We consider three copies of the system, A, B and C, submitted to the same noise 
{piiP^j- The initial condition for A is random, while it is ordered for B (erf (0) = +1), 
and C (erf (0) = —1) where i — 1, . . . , N labels the sites. A spin of of copy A is said 
to be persistent up to time t, if it experienced thermal fluctuations only, i.e. if its 



12 



history was either that of of or that of cxf. Hence 

R(t) = P(at{t') = of (*'), or aftf) = of (f ),Vf < t) (19) 

measures the fraction of persistent spins in one of the two phases ± jl4], |l^| . We also 
checked, using the determination of the interfaces proposed in ju), the equivalence 
between the definition of R{t) given by (^) and that obtained by considering a spin 
as persistent if it was not swept by an interface. 

Figure |i"o| depicts the results. The exponent 9 seems constant, with a value w 0.22, 
for all models on the low temperature side of the critical line that we considered. It is 
nevertheless hard to conclude, on the basis of numerical measurements. The decrease 
of persistence at criticality is much faster. A possible explanation is that, as for the 



voter model (see 4.1), the thickness of interfaces diverges, leading to a behaviour 
similar to (|l7|). 



5. Conclusion 



The aim of the present work was to initiate the study of coarsening and persistence 
for a class of models interpolating between the voter and Ising models, introduced 
previously in ]l~6| ]. 

In our opinion, the introduction of two temperatures T\ and T 2 , respectively 
associated to interfacial and bulk noise, gives a useful interpretation of the phenomena 
observed when varying continuously the parameters defining this class of models. The 
present analysis seems to indicate that 9 « 0.22 is the universal persistence exponent 
for all models in the low temperature phase. The situation on the critical line is 
less clear, and possibly related to the divergence of the interfaces between domains 
at criticality. In this respect it would be desirable to find a way of computing the 
equation of the critical line in the (pi,j?2) parameter plane. A more ambitious goal 
would be to compute 9 for the two dimensional zero temperature Ising model, in the 
dual reaction diffusion framework given in the appendix, using the field theoretical 
methods of @ || || 0. 

Let us finally mention extensions to d > 2 or to the Potts 
model as possible further studies. We hope to address some of those points in the 
future. 
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Appendix A. The dual approach 

In this appendix we show that the dynamics of this class of models possesses a dual 
description in terms of reaction diffusion processes for a set of random walkers moving 
backward in time. This description, beyond its intuitive value, is the starting point 
of an analysis of the temporal evolution of the system by means of diagrammatic 
expansions |p3| . 

We first review the case of the voter model [|[ [si 0, 24 . Let us define the random 



variable r = (1 + <x)/2 = 1,0, hereafter also called spin. The value of the spin at 
time t and site x can be traced back to the value of its ancestor at time t = 0, 
located on site y. The idea is to follow the line of continuity of the value of the spin, 
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backward in time, from (x,t) to (y, 0). Indeed, one first follows the line of constant r 
backward in time, from (x, t) until the last updating of site x is met. At this time, the 
spin had chosen the value of a neighbouring site. Moving to this site one again follows 
the line of constant r backward in time. Proceeding until t — 0, one thus performs a 
random walk, backward in time, from (x, t) to (y, 0) such that — Ty . 

It is easy to see that all spins at time t possess an ancestor at time t = 0, but that 
conversely not all sites at time t = are ancestors of spins at time t. Indeed, starting 
from n spins t{ , . . . , r„ , located on sites x\, .. . ,x n , at time t, the same reasoning 
leads to the consideration of n random walkers starting from sites X\, . . . ,x n at time 
t, and going backward in time. When two such walkers meet, they coalesce because 
the corresponding sites have a common ancestor from which they inherit the common 
value of their spin. 




Figure Al. An example of coalescing random walks moving backward in time, 
occurring in the voter model. The drawing is done in one dimension, for simplicity. 



Hence the determination of the values of these n spins is equivalent, for the voter 
model, to knowing the history of n independent coalescing random walkers, starting 
from sites x\, . . . , x n , and moving backward in time. An example is displayed in figure 



Al. The final stage, at time t = 0, contains k walkers (1 < k < n) at sites yi, ... ,yk- 



Since Tx} ■ ■ ■ ri* = • • • , we conclude that the sum of the probabilities of all 
these processes starting from the n points Xi,...,x n gives the correlation function 
{tx} ■ ■ ■ Tx})- One can easily generalize the reasoning to the case of spins on different 
sites X\ : . . . , x ni at different times tx, . . . , t n . 

For the voter model, the definition (Q) can be rewritten as 

P(?X = 1) = \ {T Vl + Ty 2 + Ty 3 +T y4 ). (Al) 

This means that for a two dimensional random walk the elementary process consisting 
of a jump from site x to one of its neighbouring sites y has probability \t v . Similarly 
(|ll]) leads to 

P(T X = 1) = 1 + (P2 - Pi) {T V1 + Ty 2 + Ty 3 + TyJ , (A2) 



the interpretation of which is given below. The generalization to a generic model 
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(pi,£>2) is now straightforward. By its very definition one has 

P(t x = 1) = P2T Vl Ty 2 T y3 T y4 

+Pi [T yi Ty 2 T y3 (1 - T y J + 3 perm, terms] 

+5 i T vi T V2 (1 ~ r y 3 ) (1 - T yJ + 5 perm, terms] (A3) 
+(1 - Pi) i T yi (! - T y 2 ) (! - Tys) i 1 - T VI ) + 3 perm, terms] 

+ (1 - ft) (1 - T yi ) (1 - Ty 2 ) (1 - Ty 3 ) (1 - TyJ 

This expression can be rewritten as 

P(7i = l) = l-P2 

+ (P2 - Pi) (Tj,! + T y2 + Ty 3 + TyJ 

+ (2pi - p 2 - |) (t^t^ + t Wi t ot + r yi r y4 + (A4) 
^~ T V2 T V3 + T V2 T yi ~l~ T y3 T yi) 
~ (2pi — P2 — 2) ( T yi T y2 T y3 + T yi T y2 T y4 + T y\ T y?, T yi + T y2 T y-3 T yi) 



Equation (A4) is the generalization of ( |A2| ) and can be expressed as the sum of four 
elementary processes, namely: 

• the random walker disappears (with a weight 1 — P2), 

• the random walker makes a step to a neighbouring site (with a weight P2 — Pi), 

• the random walker splits into two walkers located on different neighbouring sites 
(with a weight 2pi — P2 — 5), 

• the random walker splits into three walkers located on different neighbouring sites 
(with a weight — 2p\ + P2 + j)- 

Hence the coalescing random walks considered above for the voter model have now 
also the possibility of annihilating, or of branching into two or three new walks. Note 
that the weight P2 — Pi vanishes on the line of the majority vote model, while the 
weight 2pi — P2 — \ vanishes on the line of the noisy voter model. 

Let us underline the fact that the weights corresponding to these elementary 
processes can be negative. Therefore these drawings should be seen as diagrams, and 
not as the actual paths of real random walkers. They may indeed be used for actual 
diagrammatic computations. For example, near the voter model, one can perform a 
double series expansion of the correlation functions in the 'coupling constants' 1 — P2 
and 2p\ ~p2 — \, in order to compute the rate of the exponential decay of persistence, 
in good agreement with numerical measurements [B3|. 
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